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_ Abstract 

o 

C/3 I The Kelvin force in a layer of magnetic fluid subjected to a homogeneous 

. magnetic field and local heating is studied. The study is motivated by the 

question about the corresponding Kelvin force density [M. Liu, Phys. Rev. Lett. 
, 84, 2762 (2000)]. It is shown that the usual and the newly proposed formulation 

^ ' of the Kelvin force are entirely equivalent. It is only when approximations are 

introduced that differences arise. 
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Two prominent directions of interest can be identified among the present studies 
of phenomena occuring in magnetic fluids (MFs). These are the study of uncon- 
ventional fluid dynamics phenomena such as the "negative viscosity" effect and the 
Weissenberg effect (for a review see Q ) and the proposal of new theoretical concepts 
[^. Both directions are interwoven in the observation of a novel convective instabil- 
ity in a horizontal MF layer ||3|, ^ and in the subsequent discussion about the correct 
form of the magnetic (or Kelvin) force density |5|, |^ . 

This first discussion about the range of validity of the Kelvin force was followed 
by a second one ||^, ^ triggered by a paper announcing that a pendulum experiment 
had confirmed the invalidation of the Kelvin force in MFs The claim was not 
accomplished according to and the resulting need for a clarification entailed an 



extended paper |1C]. A clarification as in the pendulum experiment is lacking for 
the convection experiment. The aim of the present paper is to show that the usual 
and the proposed formulation of the Kelvin force are entirely equivalent. It is only 
when one introduces approximations that differences arise. 

In |3|, Q a horizontal layer of MF (stable colloidal suspension of magnetite nanopar- 
ticles dispersed in kerosene) between two glass plates is locally heated by a focused 
laser beam. It passes perpendicularly through the layer in the presence of a homo- 
geneous vertical magnetic field. The absorption of the light by the fluid generates 
a temperature gradient and subsequently a refractive index gradient. This gradient 
is optically equivalent to a diverging lens, leading to an enhancement of the beam 
divergence. As result, depending on the strength of the magnetic field different 
diffraction patterns appear Q. 

To explain the observed phenomena, the form of the magnetic force inside the 
fluid has to be known. Therefore a horizontal layer of MF is considered which is 
subjected to a homogeneous vertical magnetic induction. Since the temperature and 
the concentration of the fluid may vary, the magnetic field in the fluid is inhomoge- 
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neous and give rise to a finite Kelvin force density, fx- It can be derived from the 
Helmholtz force 



fx = /^ograd 



2 ^dp 



jj2 

fiQ—^giadx , 



(1) 



where x = cep{l+Piap) is the susceptibihty of the MF, p its density, Hint the absolute 
value of the magnetic field inside the fluid, and /xq = 47r x 10~^ H/m. Higher order 
terms in p are included in x in order to determine the Kelvin force beyond the dilute 
limit. This limit is given by /3i = 0, i.e. X = XL which is the susceptibility according 
to Langevins theory which assumes non-interacting particles. In this approximation 
XL depends linearly on the density , xl = ctp = porn^ p / {SkTnief f) , where rueff 
is the effective mass of a ferromagnetic particle with its 'attached' carrier liquid 



molecules |12], m the magnitude of the magnetic moment of the particles, T the 
temperature, and k the Boltzmann constant. The coefficient Pi of the quadratic 
term in p was determined in different microscopic models [13, 14, ^] which all 
provide the same value /3i = 1/3. In the presence of a uniform external magnetic 
induction Be^t, the internal field is given by Hint = ^ext/ipoi^ + x))- Inserting all 
expressions in Eq. ([l|), the Kelvin force follows as 

Bit Xl {1 + A [3XL (1 + (3iXl) - 1]} gradxL 



fK(x) 



XL 



(2) 



In m a variant form for the Kelvin force is proposed. By defining a different 
susceptibility x via M = (x//Uo)Bj„( with x = x/(l + x) the Helmholtz force has 
now the variant form |^ 
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grad 
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grad 



dp 



X 
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'Bint (M) is the magnetic induction (magnetization) in the fluid and Bint its absolute 
value. Due to the uniform form of the external induction and the continuity of the 
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magnetic induction across the interface, the last term is zero and the first term gives 
with the definition of x 



fv(x) = l^grad ' - ^ 



(4) 



+ 

Executing the differentiation in Eq. (^) leads exactly to the same result as in Eq. (||) . 
Therefore both formulations are indeed physically equivalent under the inclusion of 
a quadratic term in p. The equivalence is true also for higher terms in p provided 
that the susceptibility can be written as x = ctp [l + X^i^i ft(ap)*] • Thus there is 
no a priori reason to prefer Eq. (|3|) over Eq. (|l]) because both formulations lead to 
the same result as long as the definition of x is used in Eq. . This very basic fact 
independent on the relation of x on /> has to be emphasized since it recedes in the 
wake of the discussion |5|, |6| . 

The discussion about the range of validity of fx in is based on the simul- 
taneous approximation that X ^ P X ~ The concurrent correctness of both 
relations has to be checked very cautiously. Since p = xl/oi with constant a, the pro- 
portionality to the density p is equivalent with the proportionality to the Langevin 
susceptibility xl- Restricting the dependence of the susceptibilities on xl up to the 
third order, one has 

X = XL [1 + /3iXL + hxl] + 0{xi) (5) 

and 

X - x(i-x + x^ ) 

= XL[l + (/3i-l)XL + (/32 + l-2/3i)xi] +0(xi), (6) 

where the expansion (P) is valid for x 1 only. Assuming a linear dependence of 
X on XL) i-e. /5i = /92 = 0, the expansion (|6|) implies necessarily that x depends on 
higher order terms of XL- Figure [l| shows the linear behaviour of x = XL (dashed 
line) and the nonZmear behaviour of x = Xl(1— Xl+Xl) (solid line) for < xl ^ 0.5. 



A nonlinear dependence of x on XL in the region, where X ^ XL holds, is confirmed 
also by measurements (see Fig. 1(b) in 

Prom Fig. |l| it becomes evident that in the region x = XL ^ subregion xl ^ 0.06 
exists, where additionally X = XL is fulfilled. Inserting x = XL in Eqs. (|2|,Q), the 
resulting force density 

f I ^ f , ^ B'^xtXLgradxL 

fK(x) = fv(x) = -^(^^^ (7) 

is nonzero. This agreement confirms the above general statement that the usual 
and the variant form of the Kelvin force density are equivalent provided they are 
functions of x (see Eqs. (§,§)). But inserting x = XL, M = xHmt in Eq. (|l|) and 
X = XL, M = (x//^o)Bi„t in Eq. (H), respectively, one gets 

fK = i2o{Mgrad)Uint (8) 



versus 



IV 



(Mgrad)B,„t . (9) 



Where the first expression gives a nonzero force density equal to (^), it is zero in the 
second case. The reason for the difference between the nonzero result of Eq. (0) and 
the zero one of Eq. (|9|) is the following: the variant form of the Helmholtz force (^) 
is a direct function of any approximation of x whereas the correct definition of x was 
incorporated into fv(x) (see Eq. (|^) . That is the deeper reason why approximations 
cause differences if the two formulae for the Helmholtz force are used. It has to be 
noted that this discrepancy is limited to a small subregion x = XL ^ 0.06 which 
is outside the usual experimental fluids. The lowest susceptibility of commercially 
available fluids is 0.13 [|l6| . 

For X = XL > 0.06 a truncation in the expansion of x after the linear order is 
deficient (see Fig. [l|). If one inserts instead the entire term X = Xl(1 — XL + xi) in 
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Eq. (^), one obtains a nonzero force density also for fy, 

fv = -^(1 - 3xL)XLgradxL , (10) 

which is a good approximation of (0) for small xl- 

These theoretical calculations as well as the experimental measurements in 
show apparently that (i) a linear dependence of x on the density results not neces- 
sarily in a linear dependence of x on the density and (ii) nonlinear contributions of p 
are relevant for x even in the region ^ *C 1. (iii) The two formulae for the Helmholtz 
force are entirely equivalent. It is only when one introduces approximations that 
differences arise in a small subregion, x = X = XL ^ 0.06. 

References 

[1] S. Odenbach, Int. J. Mod. Phys. B 14, 1615 (2000). 

[2] M. Liu, Phys. Rev. Lett. 70, 3580 (1993); 74, 4535 (1995); 80, 2937 (1998). 

[3] T. Du and W. Luo, Appl. Phys. Lett. 72, 272 (1998). 

[4] W. Luo, T. Du, and J. Huang, Phys. Rev. Lett. 82, 4134 (1999). 

[5] M. Liu, Phys. Rev. Lett. 84, 2762 (2000). 

[6] W. Luo, T. Du, and J. Huang, Phys. Rev. Lett. 84, 2763 (2000). 

[7] A. Engel, Phys. Rev. Lett. 86, 4978 (2001). 

[8] M. Liu, Phys. Rev. Lett. 86, 4979 (2001). 

[9] S. Odenbach and M. Liu, Phys. Rev. Lett. 86, 328 (2001). 

[10] A. Engel and R. Friedrichs, On the electromagnetic force on a polarizable body, 
to appear in Am. J. Phys., ( |cond-mat/0105265D . 



6 



[11] L. D. Landau and E. M. Lifshitz, Electrodynamics of continuous media, (Perg- 
amon, Oxford, 1984), Sec. 15 and 34. 

[12] V. G. Bashtovoy, B. M. Berkovsky, and A. N. Vislovich, Introduction to thermo- 
dynamics of magnetic fluids, (Hemisphere Publishing Corp., Washington, 1988), 
Sec. 2. 

[13] L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 

[14] M. S. Wertheim, J. Chem. Phys. 55, 4291 (1971); K. I. Morozov and A. V. 
Lebedev, J. Magn. Magn. Mat. 85, 51 (1990). 

[15] Y. A. Buyevich and A. O. Ivanov, Physica A 190, 276 (1992). 

[16] Data sheets of Ferrofluidics Corporation. 

0.5 

0.4 

0.3 

?5 

0.2 
0.1 


0.1 0.2 0.3 0.4 0.5 

Xl 

Figure 1: Nonhnear dependence of the susceptibiUty x = Xl(1 ~ Xl + Xl) (sohd 
hne) on the Langevin susceptibihty xl- The dashed hne shows the hnear function 
X = XL- 
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